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Abstract 

The effect of a random phase diffuser on ffuctuations of laser light 
(scintillations) is studied. Not only spatial but also temporal phase 
variations introduced by the phase diffuser are analyzed. The explicit 
dependence of the scintillation index on finite-time phase variations is 
obtained for long propagation paths. It is shown that for large ampli- 
tudes of phase fluctuations, a finite-time effect decreases the ability of 
phase diffuser to suppress the scintillations. 



1 Introduction 

Studies of laser beams propagating through turbulent atmospheres are im- 
portant for many applications such as remote sensing, tracking, and long- 
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distance optical communications. Howerver, fully coherent laser beams are 
very sensitive to fluctuations of the atmospheric refractive index. The ini- 
tially coherent laser beam acquires some properties of Gaussian statistics in 
course of its propagation through the turbulence. As a result, the noise/signal 
ratio approaches unity for long-distance propagation. (See, for example, 
Refs.pP-[S]). This unfavourable effect limits the performance of communica- 
tion channels. To mitigate this negative effect the use of partially (spatially) 
coherent beams was proposed. The coherent laser beam can be transformed 
into a partially coherent beam by means of a phase diffuser placed near the 
exit aperture. This diffuser introduces an additional phase (randomly vary- 
ing in space and time) to the wave front of the outgoing radiation. Statistical 
characteristics of the random phase determine the initial transverse coher- 
ence length of the beam. It is shown in Refs. [7J,[8J that a considerable 
decrease in the noise/signal ratio can occur under following conditions: (i) 
the ratio of the initial transverse coherence length, Ac, to the beam radius, 
To, should be essentially smaller than unity; and (ii) the characteristic time 
of phase variations, r^, should be much smaller than the integration time, T, 
of the detector. However, only limiting cases {rd/T) — >• and {rd/T) —>■ oo 
have been considered in the literature. (See, for example, Refs. [7], [8] and 
Ref. [9], respectively). It is evident that the inequality Td « T can be 
easily satisfied by choosing a detector with very long integration time. At 
the same time, this kind of the detector cannot distinguish different signals 
within the interval T. This means that the resolution of the receiving system 
might become too low for the case of large T. On the other hand, there is a 
technical restriction on phase diffusers: up to now their characteristic times, 
Td, are not smaller than 10~^s. Besides that, in some specific cases (see, for 
example, Ref. [ID]), the spectral broadening of laser radiation due to the 
phase diffuser (Acu ~ r^^) may become unacceptably high. 

The factors mentioned above impose serious restrictions on the physical 
characteristics of phase diffusers which could be potentially useful for sup- 
pressing the intensity fluctuations. An adequate choice of diffusers may be 
facilitated if we know in detail the effect of finite-time phase variation, in- 
troduced by them, on the photon statistics. In this case, it is possible to 
control the performance of communication systems. In what follows, we will 
obtain theoretically the dependence of scintillation index on r^/T without 
any restrictions on the value of this ratio. This is the main purpose of our 
paper. Further analysis is based on the formalism developed in Ref. [8] and 
modified here to understand the case of finite-time dynamics of the phase 
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diffuser. 



2 The method of photon distribution func- 
tion in the problem of scintillations. 

The detectors of the absorbed type do not sense the instantaneous intensity 
of electromagnetic waves I{t). They sense the intensity averaged over some 
finite interval T i.e. 

I ft+T/2 

m = / mdt'. (1) 

1 Jt-T/2 

Usually, the averaging time T (the integration time of the detector) is much 
smaller than the characteristic time of the turbulence variation, ta, {T << 
ta ~ 10~^ — 10~^s). Therefore, the average value of the intensity can be ob- 
tained by further averaging of Eq. [Hover many measurements corresponding 
various realizations of the refractive-index configurations. 

The scintillation index determining the mean-square fluctuations of the 
intensity is defined by 

(2) 

where the symbol {: ... :} indicates the normal ordering of the creation and 
annihilation operators which determine the intensity, I{t). (See more details 
in Refs. [8], [11]). The brackets < ... > indicate quantum-mechanical and 
atmospheric averagings. 

The intensity / depends not only on t, but also on the spatial variable 
r. Therefore, the detected intensity is the intensity /(r, t) averaged not only 
over t as in Eq. [H but also over the detector aperture. For simplicity, we will 
restrict ourselves to calculations of the intensity correlations for coinciding 
spatial points that correspond to "small" detector aperture. This simplifica- 
tion is quite reasonable for a long-distance propagation path of the beam. 

In the case of quasimonochromatic light, we can choose /(r, t) in the form 

nr,t) = :^Ee-"^^WWWW, (3) 

^ q,k 

where h^{t) and 6q(t) are the creation and annihilation operators of photons 
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with momentum ^q. They are given in the Heisenberg representation. V = 
L^LyLz is the volume of the system. 

It follows from Eqs. [2|3] that cr^ can be obtained if one knows the average 

( : I{r,t)I{r,t') : ) = ^ E ^-^^'^''^^(6+ ^(t')V-is:(t')V| W)- 

q,k,q 

(4) 

It is a complex problem to obtain this value for arbitrary turbulence strengths 
and propagation distances. Nevertheless, the following qualitative reasoning 
can help to do this in the case of strong turbulence. We have mentioned that 
the laser light acquires the properties of Gaussian statistics in the course of 
its propagation through the turbulent atmosphere. As a result, in the limit 
of infinitely long propagation path, z, only "diagonal" terms, i.e. terms with 
(i) k = k' = or (ii) q + k/2 = q' - k' /2, q - k/2 = q' + k' /2 contribute to 
the right part of Eq. HI For large but still finite z, there exist small ranges 
of k, k' ^ in case (i) and q + k/2 q' - k72, q - k/2 ^ q' + k72 in 
case (ii) contributing into the sum in Eq. HI The presence of the mentioned 
regions is due to the two possible ways of correlating of four different waves 
(see Ref. [12]) which enter the right hand side of Eq. HI As explained in Ref. 
|13j . the characteristic sizes of regions (i) and (ii) depend on the atmospheric 
broadening of beam radii as {ARb)~^, thus decreasing with increasing z. 
In the case of long-distance propagation, (Ai?)^^ is much smaller than the 
component of photon wave-vectors perpendicular to the z axis. The last 
quantity grows with z as z^^"^. (See Ref. [Sj). For this reason, the overlapping 
of regions (i) and (ii) can be neglected. In this case Eq. H] can be rewritten 
in the convenient form: 

(:/(r,t)/(r,t'):) = ^ E E ( ^ ^-^'^Cj^WViW >< (5) 

^ |k|,|k|'<koq,q' ' 
q +— H 2 ^ 2 1 +"2" ^ 2 

where the value ko, confining summation over k, k', is chosen to be greater 
than R^^ but much smaller than the characteristic transverse wave vector of 
the photons; this is consistent with the above explanations. The two terms in 
the right-hand side correspond to the two regions of four-wave correlations. 
The quantity 

/(r,q,t) = :^ E e"^'^^:+isWV|W (6) 

^ |k|<feo ^ 
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entering the right side of Eq. O is the operator of photon density in phase 
space (the photon distribution function in r, q space) . It was used in Refs. 
[8j,[13] and fT^ for the description of photon propagation in turbulent atmo- 
spheres. By analogy, we can define the two-time distribution function 

/(r,q,t,t') = ^ E e-'^^X^.itA^.it'). (7) 

^ |k|<fco ' 

Then Eq. [5] can be rewritten in terms of the distribution functions as 

( : /(r, t)I{r, t') : ) = { : /(r, q, t)/(r, q', t') + /(r, q, t, t')f{r, q', t' , t) : ) 
q,q' 

(8) 

Let us represent /(r, q, t, t') in the form /(r, q, t, t + r). We assume that 
r < T <^ r^, as explained in the text after Eq{T] In this case the Hamiltonian 
of photons in a turbulent atmosphere can be considered to be independent 
of time. As a result, both functions defined by Eqs. [6] and [7| satisfy the same 
kinetic equation, i.e. 

(9i + Cq5r + F9q)/(r,q,t)=0 

{dt + Cq9r + F9q)/(r, q, t, t + t) = 0, (9) 

where Cq = duJc^/d({ is the photon velocity, F is a random force, caused by 
the turbulence. This force is equal to ojodnij) / Or , where ojq is the frequency 
of laser radiation. n(r) is the refractive index of the atmosphere. The general 
solution of the equation for /(r, q, t, t -|- r) can be written in the form 

/(r,q,t,t + r) = i ^ e-'^'-^*")^;^^^)^, (to) Vo)-t + ^^^^ 

|k|<fco 

where ^ 

r(to) = ^ + jj Ht'W (11) 

The functions r(t') and q(t') obey the equations of motion 

r(t')=Cq(t'),q(t') = F(r(t')) (12) 
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with the boundary conditions r(t' = t) = r, q(t' = t) = q. The instant to is 
equal to t — 2;/c, where c is the speed of hght. to is the time of the exit of 
photons from the source. This choice of to makes it possible to neglect the 
influence of the turbulence on the initial values of operators 6+ (to), b{to + r) 
(their dependence on time is as in vacuum). 

The term for /(r,q, t) can be obtained from Eq. [12] by putting r = 0. 
Substituting both distribution functions into Eq. [HI we obtain 

( : /(r, t)/(r, t') : ) = ^ E ( ^ e-^'^'-<'(*°)"*'^'v(*o) x (13) 

^ |k|,|k|'<koq,q' 

[^iM+|(^°)^(*o)-t(^o)^i(.o)+^^ + ^)V(to)-i|^(^o + r) + 

where rq(to) and rq'(to) are solutions of Eqs. [12] with the initial conditions 
r{t' = t) = r, q(t' = t) = q and r(t' = t) = r, q(t' = t) = q', respectively. 



3 A phase difFuser with finite correlation time 

The operators on the right side of Eq. [13] are related through matching 
conditions with the amplitudes of the exiting laser radiation (see Ref. [8]) 
by the relation 

V,5o = KL.Ly)-'/^ I dr^e-^^^^^<l>(r^), (14) 

where b is the operator of the laser field which is assumed to be a single-mode 
field and the subscript (_l) means perpendicular to the 2;-axis component. The 
function $ describes the profile of the laser mode, which is assumed to be 
Gaussian-type function [$ = {^)^^'^rQ^e~'^^^^^o]. ro desribes the initial radius 
of the beam. 

To account for the effect of the phase diffuser, a factor e*'^*^'"-^'*'') or e*'^'^'"-L'*o+'^) 
should be inserted into the integrand of Eq. [TU The quantity ^p{r±,t) is the 
random phase introduced by the phase diffuser. A similar consideration is 
applicable to each of four photon operators entering both terms in square 
brackets of Eq. [13] It can be easily seen that the factor 

T = p*['^('''*o)-'^('"''*o)+¥'(ri,to+T)-'/3{r'^,io+^)] ('1 c\ 
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Figure 1: Schematics of random phase variations along the x direction. Two 
curves correspond to different instants t and t + t. 



describing the effect of phase screen on the beam, enters imphcitly the inte- 
grand of Eq. [13] (the indices ± are omitted here for the sake of brevity) . There 
are integrations over variables r,r',ri,r'j^ as shown in Eq. [TH Furthermore, 
the brackets < ... >, which indicate averaging over different realizations of 
the atmosperic inhomogeneities, also indicate averaging over different states 
of the phase diffuser. As long as both types of averaging do not correlate, 
the factor f|T5l) entering Eq. [13] must be averaged over different instants, tg. 
To begin with, let us consider the simplest case of two phase correlations 



i[ip{r,to)-tf{r',to)] 



(16) 



It is evident that in the case (Iv^l) >> 1, as shown schematically in Fig. 1, 
the factor ( |T6i) is sizable if only points r and r' are close to one another. 
Therefore, the term given by Eq. [16] can be replaced by 



(17) 



where ^'^i'"'*"^ is considered to be a Gaussian random variable with the mean 

or 



square values given by ([^^^^]') = ([^^^g^]') = 2A-^ where A, is the 
correlation length of phase fluctuations. (See Fig.l). As we see, in this 
case the effect of phase fluctuations can be described by the Schell model 

m-mm-m- 

A somewhat more complex situation is for the average value of T given 
by Eq. [TJl There is an effective phase correlation not only in the case of 



7 



coincident times, but also for differing times. For Ac << tq, two different 
sets of coordinates contribute considerably to phase correlations. This can 
be described mathematically as 

(T) = |^e«[¥'(i''*o)-¥'{r',io)+¥'(ri,io+-r)-¥'(r;,to+T)]^ ^ ^gi[vp(r,to)-</j(r',io)]^ ^ ^^g) 
(^Qi[v>{ri,to+r)-<f{r'j_,to+T)]^ _|_ |/gj[</3(r,to)-</p{r; ,to+r)] ^ ^gj[v5(ri,to+r)-</3{r',io)] ^ _ 

Repeating the arguments leading to Eq. [T71 we represent the difference in 
the last term ip{r, to) — (^(r'-L, to + t) as 

d^jrto) d^{r,to)^_ (19) 

or oto 

Then, considering the random functions ^'^^^'^^ and ^"^^f"* as independent 
Gaussian variables, we obtain a simple expression for (T). It is given by 

where {[^^f) = 

As we see, the effect of the phase screen can be described by two param- 
eters, Ac and u, which characterize the spatial and temporal coherence of the 
laser beam. In the limiting case, u — > oo, the second term in Eq. [201 vanishes 
and the problem is reduced to the case considered in Refs. [Z|,[5]. In the 
opposite case, z/ — 0, both terms in Eq. [20] are important. This is shown 
in Ref. |9j. In what follows, we will see that these two limiting cases have 
physical interpretations where where » 1 (slow detector) and << 1 
(fast detector), respectively. 

There is a specific realization of the diffuser in which a random phase 
distribution moves across the beam. (This situation can be modeled by a 
rotating transparent disk with large diameter and varying thickness.) The 
phase depends here on the only variable r — vt, i.e. 

y.(r,t) =y.(r-vt), (21) 

where v is the velocity of the drift. Then we have 

(T) ^ e-Ac-'[(r-r')2+(ri-r'i)2] ^ ^_A-2[(r-r;+vr)2+(r'-ri+vr)2] _ ^22) 

Comparing Eqs. [20] and [22|, we see that the quantity, v/Xc, stands for the 
characteristic parameter describing the efficiency of the phase diffuser. The 
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criterion of "slow" detector requires {vT/X^ » 1. Qualitatively, the two 
scenarios of phase variations, given by Eqs. [2D] and [221 affect in a similar way 
the intensity fluctuations. In what follows, we consider the first of them as 
the simplest one. (This is because the spatial and temporal variables in (T), 
given by Eq. [201 are separable.) 




z(km) 



Figure 2: The scintillation index o"^ vs propagation distance z in the case of 
"slow" detector: vT oo. The parameter r^/rg indicates different initial 
coherence length. In the absence of phase diffuser (ri/ro)^ = 1 (solid line). 

is the conventional parameter describing a strength of the atmospheric 
turbulence. 

Substituting the expressions for operators given by Eq. [TH with account 
for the phase factors e^*'^'-'"-^'*^) and averaging over time as shown in Eq. [H 
we obtain 



( : m : ) 



27rr2 
LxLyV 



|k|,|k|'<koq,q' ^ 



f23) 
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where the notation < ... > after sums indicates averaging over different re- 
ahzations of the atmospheric refractive index. The parameter = TqII + 
2rQA~^)~^ describes the initial coherence length modified by the phase dif- 
fuser. Other notations are defined by following relations 

, k + k' , k + k' 

^ _ qto + Qto , k - k^ _ qto + qlo k - k^ 

~ 2 4 ' ° ~ 2 4 ' 

Further calculations follow the scheme described in Ref [8]. Fig. 2 illustrates 
the effect of the phase diffuser on scintillations in the limit of a "slow" de- 
tector (z/T oo). We can see a considerable decrease in caused by the 
phase diffuser. At the same time, the effect of the phase screen on cr^ becomes 
weaker for finite values of z/T. Moreover, comparing the two upper curves in 
Fig. 3, we see the opposite effect: slow phase variations {vT = 1) result in 
increased scintillations. There is a simple explanation for this phenomenon: 
the noise generated by the turbulence is complemented by the noise arising 
from the random phase screen. The integration time of the detector, T, is 
not sufficiently large for averaging phase variations generated by the diffuser. 

The function, \E'(z/T), has a very simple form in the two limits: (i)\E'(z/T) ^ 
1 + J)^rp , when uT >> 1 ; and (ii) ^'(z/T) «i 2, when z/T << 1. Then, in case 
(i) and for small values of the initial coherence [(rf /rp) << 1], the asymptotic 
term for the scintillation index {z — > oo) is given by 



2 _ ^1 



a'^^ + (24) 
rg V2z/T ^ ^ 

The right-hand side of Eq. [2l] differs from analogous one in Ref. [8] by the 
value that is much less than unity but, nevertheless, can be comparable 
or even greater than (ri/ro)^. 

In case (ii), the asymptotic value of cr^ is close to unity, coinciding with 
the results of Refs. [6] and [9] . This agrees with well known behavior of the 
scintillation index to approach unity for any source distribution, provided 
the response time of the recording instrument is short compared with the 
source coherence time. (See, for example, survey [1]). 
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A similar tendency can be seen in both Figs. 3 and 4: the curves with the 
smallest uT, used for numerical calculations {uT = 1), are close to the curves 
"without diffuser" in spite of the small initial coherence length [(ri/ro)^ = 
0.2; 0.05]. It can also be seen that all curves approach their asymptotic values 
very slowly. 




Figure 3: Dependence of scintillation index on the distance z for different 
values of the parameter z/T describing diffuser dynamics. The solid curve is 
calculated for (ri/ro)^ = 1 (without diffuser). Other curves are for (ri/ro)^ = 
0.2. 



4 Discussion 

It follows from our analysis that the scintillation index is very sensitive to 
the diffuser parameters, Ac and u, for long propagation paths. On the other 
hand, the characteristics of the irradience such as beam radius, R, and angle- 
of-arrival spread, A^, do not depend on the presence of the phase diffuser for 
large values oi z. To see this, the following analysis is useful. 
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T 1 1 1 1 1 — I — I — I 1 r 
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Figure 4: The same as in Fig. 3 but for smaller initial coherence length: 
{n/rof = 0.05. 



The beam radius expressed in terms of the distribution function is given 

by 

p2/,^ _ Eq/c?r^ri(/(r,z,q,t)) 

^ ^" E,/rfr^(/(rx,.,q,t))- ^''^ 

Straightforward calculations using Eq. [10] with r = (see Ref. [8]) result in 
the following explicit form: 

i?^(z) = f + ^ + 4T.^ (26) 
2 g^rf 

where T = 0.558C^/o ^^'^ is the inner radius of turbulent eddies, which 
in our previous calculations was assumed to be equal 27r x 10~^m. As we see, 
the third term does not depend on the diffuser parameters and it dominates 
when z oo. 

A similar situation holds for the angle-of-arrival spread, AO. (This phys- 
ical quantity is of great importance for the performance of communication 
systems based on frequency encoded information [10].) It is defined by the 
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distribution function as 



^^'^'"''^^ -go^E,(/(r.,.,q,t))- 

Simple calculations, which are very similar to those while obtaining R^, result 
in 

2 4z^ 
[Ae{r^ = 0, z)f = ^ + l2Tz - -^^{rf + ?>Tqlzf. (28) 

For long propagation paths, Eq. [28] reduces to 3Tz, which like i?^ does not 
depend on the diffuser parameters. 

As we see, for large distances z, the quantities and A6' do not depend 
on Ac and v. This contrasts with the case of the scintillation index. So 
pronounced differences can be explained by differences in the physical na- 
ture of these characteristics. It follows from Eq. [2] that the functional, cx^, 
is quadratic in the distribution function, /. Hence, four- wave correlations 
determine the value of scintillation index. The main effect of a phase dif- 
fuser on 0"^ is to destroy correlations between waves exited at different times. 
(See more explanations in Ref. [Ej). This is achieved at sufficiently small 
parameters Ac and z/~^. 

In contrast, and depend on two wave-correlations, both waves 
being given at the same instant. Therefore, the values of B? and do 
not depend on the rate of phase variations [v does not enter the factor (|T711 
describing the effect of phase diffuser]. Moreover, these quantities become 
independent of Ac at long propagation paths because light scattering on at- 
mospheric inhomogeneities prevails in this case. 

The plots in Figs. 3 anf 4 show that the finite-time effect is quite sizable 
even for very "slow" detectors {vT = 20). Our paper makes it possible to 
estimate the actual utility of phase diffusers in several physical regimes. 



5 Conclusion 

We have analyzed the effects of a diffuser on scintillations for the case of 
large-amplitude phase fluctuations. This specific case is very convenient for 
theoretial analysis because only two parameters are required to describe the 
effects of the diffuser. Phase fluctuations may occur independently in space as 
well as in time. Also, our formalism can be applied for the physical situation 
in which a spatially random phase distribution drifts across the beam. (See 
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Eq. [221) Our results show the importance of both parameters, Ac and i/, on 
the abihty of a phase diffuser to suppress scintillations. 
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